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In order to examine the rotational effect around neutron star in tensor-vector-scalar (TeVeS) 
theory, we consider the slowly rotating relativistic stars with a uniform angular velocity. As a result, 
we find that similar to the case in general relativity (GR), the angular momentum is proportional to 
the angular velocity. Additionally, as the value of coupling constant K becomes higher, the frame 
dragging in TeVeS becomes quite different distribution from that in GR, where we can also see the 
■ deviation even in the interior of star. While with smaller value of K, although the frame dragging 

approaches to that expected in GR, the induced vector field due to the rotation does not vanish 
and still exists. Thus, through the observations associated with relativistic object, one could be 
possible to distinguish the gravitational theory in strong field regime even in the case that the value 
of coupling constant K is quite small. 
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K ■ I. INTRODUCTION 

In the weak gravitational field such as solar system, there are many experiments and the validity of general rela- 
tivity (GR) has been shown, while the gravitational theory in the strong field regime are still unconstrained by the 
observations. However, the development of technology will be possible to observe compact objects with high accuracy 
and those observations could be used as a direct test of the gravitational theory in strong field regime [1|. In fact, 
there are attempts to test the gravitational theory by using surface atomic line redshifts or gravitational waves 
emitted from the neutron stars 0]. In these works, they suggest the possibility to distinguish the scalar-tensor theory 
y—i proposed in [J] from GR. 

As an alternative gravitational theory, the tensor- vector-scalar (TeVeS) theory has attracted considerable attention, 
which is proposed originally by Bekenstein Q as a covariant theory for modified Newtonian dynamics @ (see 0] for 
review of TeVeS). The advantage to adopt this theory is possible to explain the galaxy rotational curves and the Tully- 
Fisher law without the existence of dark matter Q. Additionally, TeVeS has successfully explained not only strong 
gravitational lensing Q but also the galaxy distribution through an evolving Universe without cold dark matter 
Q. In the strong gravitational region of TeVeS, the Schwarzschild solution was found by Giannios and the 
Reissner-Nordstrom solution was found by Sagi and Bekenstein Furthermore, Lasky et al. derived the Tolman- 
Oppenheimer-Volkoff (TOV) equations in TeVeS and they produced the static, spherically symmetric neutron star 
models in TeVeS [Hj]. 

Recently, there are some suggestions to distinguish TeVeS from GR by using some observations, i.e., with the redshift 
of the atomic spectral lines radiated from the surface of neutron star [12[ , with the Shapiro delays of gravitational 
waves and photons or neutrinos [r| , and with the emitted gravitational waves from neutron stars [3] ■ I n this article, 
we examine the different way to distinguish TeVeS from GR, i.e., the rotational effect around neutron stars. For this 
purpose, we consider the slowly rotating neutron stars with a uniform angular velocity. The analysis of slowly rotating 
neutron star in GR has originally done by Hartle [l5[, and subsequently many works have done in this field. Actually, 
taking into account the rotational effect is more natural and more important. For example, due to the rotational 
effect, the new oscillation family, i.e., so-called r mode, could be excited [l^lrlj] . So, this article would become a first 
step to consider the rotational effect in more complicated system in TeVeS. 

This article is organized as follows. In the next section, we describe the fundamental parts of TeVeS and the 
equations to produce the non-rotating relativistic stellar models in TeVeS, which are corresponding to TOV equations 
in GR. In section Mil we derive the equations to describe the slowly rotating relativistic stars. In order to discuss 
about the slowly rotating star models in TeVeS, as mentioned later, two variables are required, i.e., one is the variable 
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shown the frame dragging and the other is the vector field induced by the rotation. In this section, we also show 
the numerical results for many stellar models. At last, we make a conclusion in section ITVl In this article, we adopt 
the unit of c = G = 1, where c and G denote the speed of light and the gravitational constant, respectively, and the 
metric signature is (— , +, +, +). 

II. STELLAR MODELS IN TEVES 
A. TeVeS 

In this section, we briefly describe the fundamental parts of TeVeS, which are necessary for the calculations in this 
article (see for the details of this theory in Q). TeVeS is constructed with three dynamical gravitational fields, i.e., 
an Einstein metric a timelike 4- vector field W, and a scalar field ip, in addition to a nondynamical scalar field a. 
The vector field fulfills the normalization condition, g llv U^U v = — 1, and the physical metric is given by 

9»u = e- 2<p g^ - 2K„U V sinh(2(^), (2.1) 

along which usual matter such as electromagnetic fields propagate. Hereafter, all quantities in the physical frame 
are denoted with a tilde, and any quantity without a tilde is in the Einstein frame. The total action of TeVeS, S, 
contains contributions from the three dynamical fields and a matter contribution, i.e., the Einstein- Hilbert action, 
S g , the vector field's action, S v , the scalar's action, S s , and the matter action, S m . These four parts of action can be 
written down as 




where U[ a m = U a ,0 — I4p,a, h a/3 = g a/3 — U a U^ ', k and K are positive dimensionless parameters associated with the 
scalar and vector fields respectively, F is a dimensionless free function, £ is a constant length scale, A is a spacetime 
dependent Lagrange multiplier, / denotes the field variables collectively, and the covariant derivative denoted by | is 
taken with respect to g a j3 (see @ for details) . 

Although the free function, F, is not predicted by the theory, as mentioned later our results in this article are 
independent of this function and therefore our results are also independent from the value of t. The scalar field 
coupling, fc, has been constrained by Bekenstein [5( to be k ~ 0.03 by using planetary motions in the outer solar 
system. While, the restrictions on the vector field coupling, K, are less severe, but probably one should consider the 
range of < K < 2, because for K > 2 one can show that the stellar pressure diverges from the stellar center outward 
and one can not construct the stellar models [12j. Additionally, the Lagrange multiplier, A, can be described as a 
function of field variables from the vector equation with the normalization condition (see Eq. (j2.12j )). Furthermore, 
we should notice about the propagation speed. That is, Bekenstein showed that TeVeS allows for superluminal 
propagation of tensor, vector, and scalar perturbations when p < Q. Using perturbations of the various fields in 
the physical frame, this statement was shown and was independent of the matter content of the model. Thus, in order 
to avoid the violating causality, the scalar field should not be negative. Namely, one should construct the neutron star 
model keeping that the scalar field be everywhere greater than or equal to zero. With this condition, the cosmological 
value of scalar field, <p c , can be constrained as <p c ^ 0.001 [l2j|. 

By varying the total action, with respect to g 1 *", one can obtain the field equations for the tensor field 

+ e^, (2.6) 

where T^„ is the energy-momentum tensor in the physical frame, T a ^lA v ) = T ai lA v + T au U^ and G^ v is the Einstein 
tensor in the Einstein frame. Conservation of energy-momentum is therefore given in the physical frame as V = 0. 
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The other sources in Eq. (12.6p are given by 



--a 2 



(2.7) 



6^ =K \g af) U [aM U [fi , v] - -^ s g aP U [lM U [m g^ - \UJU V . (2.8) 

In the same fashion, by varying 5 with respect to TA^ and ip, one obtains the field equations for the vector and scalar 
fields; 

KUl a >fl }l } + XU a + 87rGer 2 WV/3<? Q 'V,7 = 8ttG (l - e" 4 *') g^U^p, (2.9) 
[^(ke 2 h^ W>u )h a ^, a ] {fj = kG [g a P + (1 + e-^) U a U } f af3 , (2.10) 

where (i(x) is a function defined by 2/i_F(/i) + pL 2 dF(p) / d[i = —2x. With this function p,, the nondynamical scalar 
field a is determined by 

kGo 2 = ii(ke 2 h af3 tp ta cp, ). (2.11) 

At last, the field equations in TeVeS are Eqs. flU} and (JZH) - IpH) . 

Now with the normalization condition for U^, from the vector equation (|2.9p one can calculate the Lagrange 
multiplier A; 



A = KU a U [a ' M ,fi + ^U a U p tp, a tp^ - 8ttG (1 - e -4 *") U a U e f a p. (2.12) 



Since it has been shown in the strong-field limit that p = 1 is an excellent approximation [10 . Ill| , we concentrate on 
the case with p = 1 This implies from Eq. (|2.1ip that a 2 = l/(kG). Moreover, since one can show that with 
fi = l the contribution of F to the field equations vanishes [H, |T(J [HI , our results are independent of this function 
and we drop it from the remaining discussion. At the end of this section, we should mention about the assumption 
that fi = 1. As described later, to determine the frame dragging effect we should impose the asymptotical flatness. 
Strictly speaking, in the asymptotic region, the assumption that p = 1 might not be good, but we set the numerical 
boundary to be r = 300M in this article and the results are independent from the position of numerical boundary if 
the numerical boundary set to be far from r = 300M. We consider that this numerical boundary might not be so far 
from star and as a first step we assume that p = 1 in the whole numerical region. 



B. Non-rotating Relativistic Stellar Models in TeVeS 

The equilibrium configurations of non-rotating relativistic stars in TeVeS have been investigated by Lasky et al. 
[lS^ . where they derived the Tolman-Oppenheimer-Volkoff (TOV) equations in TeVeS. In this subsection, we briefly 
show the TOV equations. Static, spherically symmetric metric can be expressed as 

ds 2 = g aP dx a dx p = ~e v(r) dt 2 + e c(r) dr 2 + r 2 dd 2 + r 2 sin 2 9d(f> 2 (2.13) 

where e -1 ^ = 1 — 2m(r)/r. In general, the vector field for a static, spherically symmetric spacetime can be described 
as — (U t (r), U r (r), 0, 0). But Giannios [lOj showed that in vacuum, the parameterized post-Newtonian (PPN) 
coefficients for a spherically symmetric, static spacetime with a non-zero W can violate observational restrictions. 
Thus in this article we only consider the case where W — 0, which is the same assumption in [l2|, EH- I n this case, 
the vector field can be fully determined from the normalization condition, such as — {er v l 2 , 0, 0, 0) . With this 
vector field, the physical metric is 

dS 2 = g af) dx a dx = -e u+2v dt 2 + e^ 2ip dr 2 + e^r 2 (d9 2 + sin 2 6d<j) 2 ) (2.14) 

We further assume the stellar matter content to be a perfect fluid, i.e., = (p + P)v,pUp + Pg^u, where p is the 
energy density, P is the pressure, and is a fluid four- velocity given by — e^U^. Then one can show that the full 
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system of equations reduces to 
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where a prime denotes a derivative with respect to r. Adding an equation of state (EOS), this system of equations 
can be closed. The stellar radius in physical frame, R, is determined by R = e^ v ^ Vs ^r s , where r s is the position of 
the stellar surface defined as the point where P — 0. As additional physical properties, one can introduce the total 
Arnowitt-Deser-Misner (ADM) mass, Madm, and the scalar mass, M v , which are defined by 
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r 2 (p + 3P) e^+^-^dr, 



(2.20) 
(2.21) 



where m^ and ip c are the mass function evaluated at radial infinity and the cosmological value of the scalar field, 
respectively. 

In particular, in this article, we adopt the same EOS as in 0, [H, E3|i which are polytropic ones derived by fitting 
functions to tabulated data of realistic EOS known as EOS A (soft EOS) and EOS II (intermediate EOS). In practice, 
the adopted EOS can be described as 
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(2.22) 

(2.23) 

(2.24) 
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where T = 2.46 and K = 0.00936 for EOS A and T — 2.34 and K = 0.0195 for EOS II. Additionally, we adopt the 
values of parameters k and ip c as k = 0.03 and ip c = 0.003, which are same choices as in while the value of K is 
considered in the range of < K < 2 (see [12| for discussion about the range of value of K). 



III. SLOWLY ROTATING RELATIVISTIC STARS IN TEVES 



A. Rotational Dragging 

In this article, we consider a slowly rotating relativistic stellar models jwith a uniform angular velocity ft, where 
we assume to keep only the linear effects in the angular velocity (see in 15] for the discussion about slowly rotating 
relativistic stars in GR). Here we put the rotational axis to be 9 = 0. In this case, the star is still spherical because 
the deformation due to the rotation is of the order J7 2 , and the metric in physical frame is given by 

ds 2 = -e u+2v dt 2 + e^^dr 2 + r^e' 2 * (d6 2 + sin 2 6d(/) 2 ) - 2ujr 2 e~ 2ip sin 2 9dtd<p, (3.1) 

where the last term in the right hand side corresponds to the rotational effect and in general lo ~ 0(fl) can be 
expressed as 

W (r,9) = -^P/, (3.2) 
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where P( — Pi(cos9) is the Legendre polynomial of order I, While, the rotational effects in the other components 
of metric become of the order SI 2 , because those should behave in the same way under a reversal in the direction of 
rotation as under a reversal in the direction of time. Similarly, the deviations of pressure, density, and scalar field 
due to the rotation should be of the order (l 2 . Hereafter, in order to distinguish the effects due to the rotation, the 
deviation from the background properties would be expressed by using variables with S. Then the fluid velocity in 
physical frame of the order SI can be described as 



Su" = (0,0,0,Qu*), 



(3.3) 



where W* is i-component of four velocity in the case without rotation, i.e., u 1 = e v v l 2 ^ and non-zero component of 
the energy-momentum tensor of the order S7 is only JT^, which are given by 

Sf t! p = -fp + pj r 2 e - 2ip (ajsm6d e Pi + nsm 2 6^ + Pujr 2 e~ 2ip sm6d 9 Pe. (3.4) 

On the other hand, in Einstein frame the deviation of metric from the spherical symmetry is determined by 

5g^ = e 2 ^%„ + (e 4 ^ - 1) {8UJA V + U^SU,) . (3.5) 

By using the normalization condition for the vector field, one can show that SU l — SUt = 0. So the deviation of vector 
field from the spherical symmetry can be described as 



SUP = [ 0, W(r)P e , V(r)d e P e , ^-^dePi 

sin 



(3.6) 



However, the variables of W and V should be of the order Sr because these components are corresponding to the 
polar parity. Then the non-zero component of metric in Einstein frame of the order SI is only 5g t( p, which is given as 



Sg t4> = r 2 \e~ 4v u + e v ' 2 (e' 4 * - l) V 
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Then from the (t, <fi) component of the Einstein equation (|2.6[) . one can get the following equation 
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On the other hand, the additional equation to express the relation between uj and V can be obtained from the 
component of the vector field equation (|2.9p . which is 
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Especially, in this article we adopt the case for I = 1 to reproduce the GR limit of TeVeS [Hj]. In fact, the other value 
of £ can not be chosen to satisfy the regularity condition at the stellar center and the asymptotic behavior far from 
the star. Then, by combining the above equations (|3.8j) and (|3.9j) . one can obtain the following equations 
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(3.11) 



Note that the GR limit of Eq. (|3.11[) agrees with the well-known equation describing the frame dragging in GR. With 
these two equations (|3. 101) and p. lip and with the appropriate boundary conditions, the distributions of ui(r) and 
V(r) can be determined. Additionally, with an asymptotic flatness, one can show that u> and V are decreasing as 1/r 3 
far from the central object. 



B. Numerical Results 



In order to determine the distributions of ui(r) and V(r) with the fixed value of K for the adopted stellar model, 
we impose two boundary conditions, i.e., the regularity condition at the stellar center and the asymptotic flatness far 
from the star. In practice, at stellar center we set that w(0) = ujq, w'(0) = 0, V(0) = Vo, and V'(0) = 0, where ujq 
and Vo are some constants. Then we find the correct values of ujq and Vo in such a way that the solutions of u}(r) 
and V(r) should satisfy the asymptotic behavior as mentioned in the previous section, i.e., with the trial values of 
Wo and Vo the above ordinary differential equations can be integrated outward and we search the correct values by 
changing the trial values iteratively, where for integration we use the 4th order Rungc-Kutta method. In this way, 
we can determine numerically the distributions of uj{r) and V(r) and those behaviors far from the central star can be 
described as 

«(r) = ^+e>(£), (3.12) 

VW = ^ + o(^). (3-13) 

It should be noticed that in GR the distribution of ui(r) outside the star can be described analytically as w(r) — 2 J/r 3 , 
where the constant J is corresponding to the total angular momentum of the star [l5j . Additionally, it is well known 
that the angular momentum in GR is linearly related to the angular velocity for slow rotation as J = where the 
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constant of proportionality / corresponds to the relativistic generalization of momentum of inertia for slowly rotating 
systems [HJ]. While, in the case of TeVeS, we can see the same feature as the case of GR, i.e., as shown in Figs. [T]and 
EJ the values of J/Q and V c /& are independent of the value of angular velocity tl if only one chooses the non-rotating 
stellar models. 



GR 

- s: = o.i 

- K = 0.2 

- a: = 0.3 - 



K = 0.45 

™0 4 8 12 16 20 

Q [kHz] 

FIG. 1: Dependence of J/Q on the angular velocity f2 for the stellar model with EOS A and Madm = 1.4Mq, where the solid 
line corresponds to the case of GR while the broken lines are results in TeVeS with different values of K. 
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""O 4 8 12 16 20 
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FIG. 2: Dependence of V c /f2 on the angular velocity Q for the stellar model with EOS A and Madm = 1.4Mq. 

Fig. [3]shows the distributions of cj(r) with different values of K for the stellar model with EOS A, Madm = 1.4M©, 
and Q = 1 kHz. At a glance, we can recognize the difference of distributions between in GR and in TeVeS. Especially, 
with higher value of K those differences become obvious even in the interior of the star. As mentioned before, the 
behavior of w(r) far from star is proportional to r~ 3 , which can be seen in the right panel of Fig. [3] On the other 
hand, Fig. [4] shows the distributions of V(r) with different values of K for the same stellar model as in Fig. [3] In the 
right panel of Fig. [4] one can see the behavior of V(r) far from star. From Figs. [3] and 01 we can find the important 
point that as the value of K becomes smaller, the distribution of ui(r) in TeVeS approaches to that in GR while the 
induced vector field V(r) does not vanish. Namely, with smaller value of K the physical metric in TeVeS is almost 
same as that in GR, but there still exists the additional vector fields induced by the stellar rotation. So if one will 
consider the stellar oscillations and/or the emitted gravitational waves in TeVeS, this induced vector field could play 
a role as source term in the linearized equations which is corresponding to the linearized right hand side in Eq. (|2.6p , 
and as a result the deviation in frequencies depending on the gravitational theory can be seen even if the physical 
metric is not so different from each other. Thus such kinds of observation could tell us the gravitational theory in the 
strong field regime even if the value of K is quite small. 

Furthermore, the dependences of J and V c on the value of K are shown in Figs. [5] and [HI where the stellar masses 
are fixed to be Madm = 1.4M@ and the angular velocity are adopted that Ct = 1 kHz for EOS A and EOS II. In Fig. [5] 
the horizontal broken lines denote the values of J in GR for each EOS. From this figure, we can see that depending on 
the value of K, the angular momentum J could become 42.7% smaller for EOS A and 65.0% smaller for EOS II than 
those expected in GR. Even if this difference cannot be directly observed, one could find a kind of evidence depending 
on the gravitational theory in the strong field regime by observing the phenomena around rotating compact object. 



80 



CJ 



60 



CS 



-40 



-50 



8 




10-3 
lO- 4 
S 10-5 

3 

10-6 
lO"? 



1 



r [km] 



10 

r [km] 



100 



FIG. 3: Distributions of co{r) with different values of K for the stellar model with EOS A, Madm = 1.4Mq, and Cl — 1 kHz. 
In order to compare the results, the solid line shows the distribution in GR. 
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FIG. 4: Distributions of V(r) with different values of K for the stellar model with EOS A, Madm 
where the different lines correspond to the results with values of K shown in Fig. [3] 



1.4M , and fi = f kHz, 



IV. CONCLUSION 

In this article, in order to examine the rotational effect around the neutron star in the tensor-vector-scalar (TeVeS) 
theory, we consider the slowly rotating relativistic stars with a uniform angular velocity Q. To deal with this problem 
in TeVeS, one has to take into account not only the usual frame dragging u> but also the induced vector field V, which 
is corresponding to the <fi component. The equations for u(r) and V(r) are derived from the Einstein and vector 
field equations, and then the distributions of those variables are determined numerically with appropriate boundary 
conditions. 

As a result, we find that, similar to the case in GR, the value of J/Cl in TeVeS is constant if the stellar mass 
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FIG. 5: Dependence of J on K for the stellar models with A^adm = 1.4M and Q = l kHz for EOS A and EOS II. The 
horizontal broken lines are corresponding to the values for each EOS in GR. 
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FIG. 6: Dependence of V c on K for the stellar models with Madm = 1.4M© and 9, = 1 kHz for EOS A and EOS II. 



and value of K are fixed, where J corresponds to the angular momentum. Additionally, with higher value of K, the 
distribution of uj(t) deviates obviously from that in GR due to the existence of induced vector field and this deviation 
can be seen even in the interior region of star. On the other hand with smaller value of K, although uj approaches to 
the that in GR, the induced vector field dose not vanish. That is, even if the stellar properties in TeVeS with small 
value of K, such as the mass and radius, are almost same as those in GR, the induced vector field in TeVeS could still 
exist. This is a crucial difference depending on the gravitational theory. Thus through the observable phenomena, 
such as the stellar oscillation and emitted gravitational waves, it is possible to distinguish the gravitational theory in 
the strong field regime. Furthermore, through the stellar magnetic effect, one might be see another effect coupling to 
the induced vector field. For example, some quasi-periodic oscillations are observed recently in the giant flares and 
these are believed to be related to the oscillations of strong magnetized neutron stars [18|. Considering the stellar 
magnetic fields, one might be possible to obtain the further constraint in the theory. 



Acknowledgments 

We thank Demetrios Papadopoulos, Kostas D. Kokkotas, Miltos Vavoulidis, and Erich Gaertig for valuable com- 
ments, and also the referees for fruitful comments. This work was supported via the Transregio 7 "Gravitational Wave 
Astronomy" financed by the Deutsche Forschungsgemeinschaft DFG (German Research Foundation). 



[I] D. Psaltis, Living Rev. Relativity 11, lrr-2009-9 (2009). 

[2] S. DeDeo and D. Psaltis, Phys. Rev. Lett. 90, 141 101 (2003). 

[3] H. Sotani and K.D. Kokkotas, Phys. Rev. D 70, 084026 (2004); 71, 124038 (2005). 

[4] T. Damour and G. Esposito-Farese, Classical Quantum Gravity 9, 2093 (1992). 

[5] J.D. Bekenstein, Phys. Rev. D 70, 083509 (2004). 

[6] M. Milgrom, Astrophys. J. 270, 365 (1983). 

[7] C. Skordis, Class. Q. Grav. 26, 143001 (2009). 

[8] D.M. Chen and H.S. Zhao, Astrophys. J. 650, L9 (2006). 

[9] S. Dodelson and M. Liguori, Phys. Rev. Lett. 97, 231301 (2006). 
[10] D. Giannios, Phys. Rev. D 71, 103511 (2005). 
[11] E. Sagi and D. Bekenstein, Phys. Rev. D 77, 024010 (2008). 
[12] P.D. Lasky, H. Sotani, and D. Giannios, Phys. Rev. D 78, 104019 (2008). 
[13] S. Desai, E.O. Kahya, and R.P Woodard, Phys. Rev. D 77, 124041 (2008). 
[14] H. Sotani, Phys. Rev. D 79, 064033 (2009); 80 064035 (2009). 
[15] J.B.Hartle, Astrophys. J. 150, 1005 (1967). 

[16] M. Vavoulidis, K.D. Kokkotas, and A. Stavridis, Mon. Not. R. Astron Soc. 384, 1711 (2008). 
[17] E. Gaertig and K.D. Kokkotas, Phys. Rev. D 80, 064026 (2009). 

[18] H. Sotani, K.D. Kokkotas, and N. Stergioulas, Mon. Not. R. Astron Soc. 375, 261 (2007); 385, L5 (2008); H. Sotani, A. 

Colaiuda, and K.D. Kokkotas, ibid. 385, 2161 (2008); H. Sotani and K.D. Kokkotas, ibid. 395, 1163 (2009). 
[19] When one will consider the phenomena on cosmological scales, the assumption that fi = 1 might not be a good choice [5J . 



